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ABSTRACT 

We numerically construct slim, global, vertically integrated models of optically 
thin, transonic accretion discs around black holes, assuming a regularity condition at 
the sonic radius and boundary conditions at the outer radius of the disc and near 
the black hole. In agreement with several previous studies, we find two branches of 
shock-free solutions, in which the cooling is dominated either by advection, or by 
local radiation. We also confirm that the part of the accretion flow where advection 
dominates is in some circumstances limited in size: it does not extend beyond a certain 
outer limiting radius. New results found in our paper concern the location of the 
limiting radius and properties of the flow near to it. In particular, we find that beyond 
the limiting radius, the advective dominated solutions match on to Shapiro, Lightman 
& Eardley (SLE) discs through a smooth transition region. Therefore, the full global 
solutions are shock-free and unlimited in size. There is no need for postulating an extra 
physical effect (e.g. evaporation) for triggering the ADAF-SLE transition. It occurs 
due to standard accretion processes described by the classic slim disc equations. 
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1 INTRODUCTION 

Most of our knowledge about accretion disc structure comes 
from a very successful, but simplified, approach introduced 
by Shakura (1972), Pringle & Rees (1972), Shakura & Sun- 
yaev (1973) and others. In this approach, only the radial 
structure of the disc is studied in a detailed way, and all of 
the disc properties are taken to be independent of time and 
azimuthal angle and are averaged in the vertical direction. 
Let H = H{r) be the vertical extension of the disc at a par- 
ticular radial location r. In the standard Shakura-Sunyaev 
'thin' model, terms of the order [H/r)^ and higher are ne- 
glected. This allows the problem to be reduced to a set of 
linear algebraic equations which it is easy to solve: the stan- 
dard model gives explicitly all of physical characteristics of 
the disc. 

However, the (H/r)^ terms, which are neglected in the 
standard thin disc model, describe efi'ects that in many as- 
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trophysical applications are very relevant: the radial pres- 
sure gradient, non-Keplerian angular momentum distribu- 
tion, and cooling by radial advection of heat. These terms 
are all retained in the 'slim' accretion discs models intro- 
duced by Abramowicz et al. (1988, hereafter ACLS88) and 
based on a set of non-linear radial differential equations de- 
rived by Paczyhski & Bisnovatyi-Kogan (1981). Mathemat- 
ically, integrating the slim disc equations leads to a non- 
standard eigenvalue problem which is surprisingly difficult 
to solve. The original ACLS88 slim disc models, and slim 
models studied later by many other authors, have been con- 
structed with a help of a very accurate pseudo-Newtonian 
model of the black hole gravitational field (Paczyriski & Wi- 
ita 1980). Lasota (1994) derived the full set of slim accre- 
tion disc equations in Kerr geometry, i.e. in the gravitational 
field of a rotating black hole. Abramowicz et al. (1996, here- 
after ACGL96), and Abramowicz, Lanza & Percival (1997) 
later made some improvements and corrections to Lasota's 
equations. Recently, Beloborodov, Abramowicz & Novikov 
(1997) rederived the relativistic slim accretion disc equations 
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again, taking into account the relativistic effect of the heat 
inertia that is important for discs around rapidly rotating 
black holes, but was neglected by previous investigators. 

Initially, all of the slim models were constructed under 
the assumption of the disc being very optically thick but 
a few years ago it was realized by Narayan & Yi (1994), 
Abramowicz et al. (1995), and Narayan & Yi (1995) that a 
new class of very hot, optically thin, advective discs exists. 
The new class are called ADAFs and are very promising 
for explaining properties of X-ray transients, low luminosity 
active galactic nuclei, the Galactic Centre, and other high 
energy objects. For this reason they are now being studied 
intensively by many researchers (see Narayan 1997, for a 
review and references). 

The first ADAF models were local (in the sense that 
no boundary conditions were applied), and were constructed 
with many simplifying assumptions. Only very recently have 
several authors constructed more accurate, slim global mod- 
els. Chen, Abramowicz & Lasota (1997), Narayan, Kato 
& Honma (1997), and Nakamura et al. (1997) constructed 
global transonic models for ADAFs using the Paczyhski- 
Witta model for the black hole gravity. ACGL96 and 
Jaroszyhski & Kurpiewski (1997) constructed global tran- 
sonic ADAF models in Kerr geometry, using Lasota's equa- 
tions. 

In this paper we construct some new black hole models 
of slim discs, using Lasota's equations in the version given by 
ACGL96 and Abramowicz et al. (1997a), where the influence 
of the heat inertia (discussed in Beloborodov et al. 1997) was 
neglected. Equations, boundary conditions and the numeri- 
cal procedure are presented in Section 2. In Section 3 we give 
the numerical results for adiabatic and non-adiabatic mod- 
els. Section 4 contains the final discussion and conclusions. 
The new type of global transonic models found in our paper 
consist of an ADAF in the iimer part and an SLE (Shapiro, 
Lightman & Eardley 1976) type of flow in the outer part. 
The SLE discs are known to be violently unstable (Piran 
1978), and therefore our solutions with the ADAF-SLE tran- 
sitions are not astrophysically realistic models of stationary 
accretion flows. However, they provide a non-trivial example 
of the ADAF- (thin disc) transition and for this reason they 
could be of interest. 
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where the pressure p is at the equatorial plane, and 
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The definitions of the symbols introduced above are given 
later in this Section. The last term in equation (2.2) differs 
from the one used by ACGL96. 
The equation of angular momentum conservation. 
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is the specific angular momentum. 

The equation of energy conservation in the one temperature 
(Te — Ti = T) approximation. 
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where 7g is the adiabatic index of the gas. 
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is the surface heat generation rate, and is the radiative 
cooling flux. For the case of proton-electron bremsstrahlung 
cooling, the emissivity per unit volume takes the form 
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The equation (2.6) is written for a gas-pressure dominated 
medium with the equation of state 

p = -pT, (2.9) 

where TZ is the gas constant, p is the mean molecular weight 
of gas (p = 1/2 for a hydrogen plasma). 
We take the equation of vertical balance in the form given 
by Abramowicz, Lanza & Percival (1997), 
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2 METHOD OF SOLUTION 
2.1 Equations 

The equations for Keplerian discs in the Kerr metric were 
derived by Novikov & Thorne (1973). We follow the nota- 
tion and equations in ACGL96 for non-Keplerian slim discs, 
except where explicitly indicated. 
The mass conservation equation reads, 



M = -27rA 



1/2, 



V 



(2.1) 



where E — 2Hp is the surface density, H is the half-thickness 
of the disc, p is the rest mass density at the equatorial plane 
and V is the radial velocity, defined later by the equation 
(2.16). 

The equation of radial momentum conservation takes the 
form. 



(2.10) 
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In equations (2.1) - (2.11), we have used the following 
definitions. 
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where M is the mass and a is the total specific angular 
momentum of the Kerr black hole. 
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are the angular velocities with respect to a stationary ob- 
server and to a local inertial observer, and u' are compo- 
nents of the four- velocity u' of matter. The angular frequen- 
cies of the corotating (+) and counterrotating (— ) Keplerian 
orbits are 



Ml/2 



Ml/2 



(2.14) 



The relation between the Boyer-Lindquist and physical ve- 
locity components in the azimuthal direction is 

i;*"^^ = m. (2.15) 

The radial velocity V is defined by the formula 

-^=u^gU\ (2.16) 

where Qrr is the metric tensor component. The velocity V is 
the radial velocity of the fluid as measured by an observer 
at fixed r who corotates with the fluid. The Lorentz gamma 
factor can be written as 
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We use the standard assumption for the viscosity coefficient, 

u = aCsH, (2.18) 

where a is a parameter and Cs ~ {p/ pY^"^ is the isothermal 
sound speed. 



2.2 Conditions at the sonic radius 

The differential equations (2.2) and (2.6) can be rewritten 
in a form where just two unknown functions V{r), and T{r) 
are present. 
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The regularity condition at the critical point Ts of equations 
(2.19) and (2.20) requires that 
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where A'^ and D are given by the expressions (2.21) and 
(2.22). In the nonrelativistic limit, the condition D = re- 
duces to V'^ — 27gCs/(7g + 1). We note that in our version of 
the vertically integrated equations, the velocity at the criti- 
cal point does not coincide with the adiabatic sound speed, 
as it does in the case of the spherical accretion (Bondi 1954). 
Nevertheless, following tradition, we use the name the 'sonic 
radius' for the critical point r^, and define the Mach number 
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2.3 Boundary conditions 

The set of equations to be solved consists of two first or- 
der differential equations (2.2) and (2.6), and one second 
order equation (2.4). To find a solution with given M, a 
and M, four integration constants (free parameters) are re- 
quired, but the regularity condition (2.23) at the sonic ra- 
dius reduces the number of free parameters in the problem 
to three. The choice of the three free parameters should be 
determined by fixing the boundary conditions. This may be 
done in many ways, but only certain particular choices will 
be consistent with stable numerical algorithms. 

The standard approach is to fix boundary conditions at 
a large distance from the central accreting black hole (e.g. 
ACLS88, Chen et al. 1997) but a modification of this was 
used by Narayan et al. (1997) who specified = fi/f and 
Cs — 10~"^nr at the outer boundary and the no-torque con- 
dition (Kl/dr = at the inner boundary located between 
the black hole horizon and the sonic radius. It is well known 
that when the standard (or modified standard) boundary 
conditions are used, the mathematical problem is an eigen- 
value problem that is numerically very time consuming and 
difficult to solve. 

To avoid this difficulty, Chakrabarti and his collabo- 
rators (see Chakrabarti 1996 for references) introduced a 
very clever mathematical trick. They assumed a different 
set of boundary conditions, which are applied not only at 
the outer boundary, but also at the sonic radius and near to 
the black hole horizon. In this way, the most difficult part of 
the problem - finding the eigenvalue - is trivially solved: the 
eigenvalue may be assumed to have any particular value (in 
practice from some continuous range) and the outer bound- 
ary conditions automatically follow from it. In a way, this is 
putting the cart before the horse: one tries to match a prob- 
lem to an already known solution. While this is obviously 
a perfectly acceptable mathematical procedure, there is no 
guarantee that solutions constructed in this way will always 
be astrophysically acceptable. One should worry about the 
appearance of unphysical discontinuities which prevent so- 
lutions extending smoothly out to very large radii. Troubled 
solutions such as these may correspond to incorrect choices 
of eigenvalues in the sense that no astrophysically acceptable 
outer boundary conditions are consistent with them. 

In order to see how qualitative differences in the treat- 
ment of the boundary conditions inffuence numerical models 
of global transonic solutions describing slim accretion discs, 
we have adopted in this paper a procedure that is very simi- 
lar to the one used by Chakrabarti and collaborators. Specif- 
ically, we take the position of the sonic radius Ts as one of 
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the free parameters. The second free parameter is fixed by 
the condition that the angular momentum tends asymptot- 
ically to a constant value near to the black hole horizon r^. 
We assume that at the inner boundary ri„ of the disc, the 
radial change of angular momentum is already negligible, 

^ =0. (2.25) 
dr 

Strictly speaking, in (2.25) one should set ri„ = Vh- In prac- 
tice, however, in numerical calculations one may choose the 
location of the inner boundary to be anywhere between the 
horizon and the sonic point, and the results change only 
slightly. We use Tin = {ru + rs)/2 in our numerical calcula- 
tions. 

The third parameter is specified at the outer boundary Tout, 
by fixing the angular velocity, 

n(ro„t) = ttout, (2.26) 
as a fraction of the Keplerian angular velocity Six (rout). 



2.4 Numerical procedure 

The set of nonlinear differential equations (2.4), (2.19) and 
(2.20) are solved numerically on a fixed numerical grid by 
means of a standard relaxation technique. This determines 
the solution by starting from an initial estimate which is 
then improved iteratively. The original differential equations 
are replaced by approximate finite-difference equations. We 
used two overlapping grids, one of which, {vi}, i = 1, ■ ■ ■ , N , 
is used for the second-order equation (2.4) while the other, 
{^1+1/2}) * = l,''',-^ — 1, is used for equations (2.19) 
and (2.20). Variables Qi, i = 1,---,N, and Vi+1/2, 7i+i/2, 
i = 1, • • • , A*' — 1, are defined on grids {r^} and {ri+1/2}, re- 
spectively. It is assumed that the position of the sonic point 
ra coincides with one of the points ^^+1/2 at i = A^s. 

The problem involves SA'^ — 6 coupled finite-difference 
(algebraical) equations which depend on 3N — 2 variables. 
This set of equations can be uniquely solved if we exclude 
two variables fli and Sljv using the conditions (2.25) and 



(2.26), respectively, and another two variables Vj 



and 



Tjv^+1/2 using the conditions (2.23) at the sonic radius. 
Then, the system of 3N — 6 algebraical equations 

Fk{xi,X2,--- ,X3N-6) = fc = l,2,---,3iV-6, (2.27) 

where 



Xk = (^^2 



, ^JV-l, V3/2, • • • , yNs-1/2, 



Vns+3/2, • • ■ , Vm-1/2, T3/2, • ■ • , Tff^_i/2, 



T, 



JVs+3/2, ■ ■ 



N-1/2 



is solved by means of a multidimensional Newton iteration 
method (see, for example. Press et al. 1992). In matrix no- 
tation, the procedure of solution for equations F(x) — can 
be represented as. 



where 

5x = -J ■ F. 



< e sC 1, 



(2.28) 



(2.29) 



Here J is the Jacobian matrix, Jij = dFi/dxj. We take a 
value of the parameter e less than 1 to obtain convergence if 
the initial guess is not sufficiently close to the solution but 
then, when x becomes close to the solution, we use the full 



Newton step, e = 1, to obtain the fastest (quadratic) conver- 
gence. Finally, we check the degree to which both functions 
and variables have converged. 



3 RESULTS 

We have constructed two families of models; a family of adi- 
abatic models in which radiative cooling was neglected and a 
family of bremsstrahlung models with bremsstrahlung cool- 
ing. In both cases we assumed a Schwarzschild black hole, 
with a = 0, and took the disc matter to consist of hydrogen 
plasma with /i = 1/2. We used two values for the adiabatic 
index in the models, 7^ = 5/3, 4/3, the first corresponding 
to gas pressure domination while the second could represent 
the case where isotropically tangled magnetic field provides 
a major contribution to the total pressure. It is often argued 
that in realistic astrophysical accretion flows, the gas and 
magnetic pressures are roughly the same. Thus, one may as- 
sume that the effective adiabatic index has a value between 
4/3 and 5/3. 



3.1 Adiabatic models 

We take = in equation (2.6) in order to study the 
pure effect of advection dominated cooling. For a given value 
of there are two main parameters of models: a and Tg. 
Solutions depend only very weakly on other parameters (see 
Section 2.3). In particular, the problem is independent of 
the mass M after a radial re-scaling with rg = 2GM/c? . 

We now describe models with adiabatic index 7g = 5/3. 
Two types of numerical solution have been found in this case. 
The first type, which are limited in size, can be constructed 
if Ts > {rs)crit, where {rs)crit depends on a. They cannot 
extend in the radial direction beyond a radius at which 
the thickness H goes to zero; it is impossible to construct 
models with rout larger than this. The value of r^t^ is a 
decreasing function of Vs- Our numerical procedure can find 
only an approximate value for rJJS^, because dH /dr becomes 
singular [dH/dr —> —00) when H —> 0, but one can set 
the outer boundary condition at rout < r^ut^ and obtain 
solutions which are just part of the solution with Tout = 
r'^ut^ ■ All of these models have the ratio H/r less than one 
everywhere, and so the condition for the slim disc approach 
to be self-consistent (Hjr ^ 1) is satisfied. Models with 
rs ~ {rs)crit have the maximal value of H/r approximately 
equal to one. Models of this type can experience a transition 
at radius « rj^"^ to a thin disc of the type allowed by the 
energy equation. In Section 3.2, an example of the transition 
of an ADAF to an SLE-type disc is given. 

The second type of solution exists when < {rs)crit- 
Models of this type are not size-limited in the radial direc- 
tion; indeed, their properties depend strongly on the value 
of rout - In particular, the maximum value of the ratio H/r is 
an increasing function of rout- Typically, these models have 
H/r ^ 1 for the whole range of radii, ri„ < r < rout (mod- 
els with H/r < 1 everywhere exist only for a small value of 
rout)- This type of solution is not self-consistent because it 
does not satisfy the slim disc condition for large r and so 
we will not discuss it further. In the following, we will only 
discuss models of the first type, which have limited radial 
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Table 1 . Parameters of the adiabatic accretion disc models (adi- 
abatic index ■jg = 5/3). The radial distances Ta, {rs)crit, ^^^t^ 
are given in units of the gravitational radius rg . 



Model 


a 


rs 


{^s )crit 


^max 
out 


Al 


0.1 


2.485 


2.485 


24.0 


A2 




2.5 




16.8 


A3 




2.55 




11.1 


Bl 


0.01 


2.165 


2.165 


85.7 


B2 




2.175 




53.3 


B3 




2.2 




29.0 


Bi 




2.25 




17.1 


B5 




2.3 




13.0 


CI 


0.001 


2.09 


2.09 


148.0 


C2 




2.1 




73.0 


C3 




2.125 




37.4 


C4. 




2.15 




25.5 


C5 




2.2 




15.7 


Dl 




2.05 


2.05 


245.0 


D2 




2.1 




35.6 


D3 




2.15 




18.6 




10 - 



2 2.2 2.4 2.6 

Figure 1. The dependence of maximal radial extension r"^^^ of 
the adiabatic accretion disc models (adiabatic index 7g = 5/3) 
on the position of the transonic radius Ts ■ The data are shown for 
the viscosity parameter a = 10" 10""', 0.01 and 0.1 (triangles, 
squares, circles and pentagons, respectively). In the region above 
the dashed line, slim disc solutions do not exist. 

Table 1 summerizes the basic parameters of the com- 
puted models with a = 10~^, 10~^, 0.01 and 0.1. In these 
models Tout was set approximately equal to r^ut°. Fig. 1 
shows the dependence of r^"'^ on the position of the sonic 
radius for the solutions listed in Table 1, with the dashed line 
showing roughly the boundary of the region in the (rs, r'^t') 
space where solutions of the first type exist. More extended 
discs correspond to smaller a and smaller r^. Fig. 2 shows 
the disc shapes H{r) for the largest discs (with rs = (rs)crit) 



10 - 



a = 0.1 




5 10 15 20 25 




50 100 150 200 250 

Figure 2. The shapes of the adiabatic accretion discs (7g = 5/3) 
for four values of the viscosity parameter a = 0.1, 0.01, 10~^ and 
10~5. The models presented have the maximal radial extensions 
for the given a, and correspond to the models Al, Bl, CI and 
Dl, listed in Table 1. 



for the values of a listed above (models Al, Bl, CI and Dl). 
The shapes of these discs are very similar: the maximum 
vertical thickness occurs at around CGSr^"^ and the outer 
edge is very sharp. Figs 3, 4, and 5 show the variation with 
r of the angular momentum £, the radial velocity V and the 
equatorial pressure p, for solutions with different values of 
a and rs (corresponding to the models listed in Table 1). 

At the outer edge, all of the models have low values 
of V and Cs, and the distribution of angular momentum 
has a super-Keplerian part. The super-Keplerian rotation 
corresponds to a region where the pressure increases with 
radius (see Fig. 5). The pressure increases due to the rapid 
decrease of the disc thickness H at the outer boundary and 
the corresponding increase in the matter density. The super- 
Keplerian rotation compensates the gravitational force and 
the gradient in the gas pressure, which are both directed 
inwards in this region. The presence of super-Keplerian ro- 
tation does not depend on the particular value taken for 
^out- We tried a variety of values, both Keplerian and non- 
Keplerian, but always found the presence of super-Keplerian 
rotation. Note, that the super-Keplerian rotation often oc- 
curs in a region where a sharp transition in qualitative 
properties of the accretion flow takes place; the two sub- 
Keplerian flows on either side of the super-Keplerian part are 
very different. For example, super-Keplerian rotation occurs 
at the sharp, cusp-like inner edge of thick discs (Abramow- 
icz, Calvani & Nobili 1980). In this case, the super-Keplerian 
part is sandwiched between supersonic free-fall with negligi- 
ble dissipation, and centrifugally supported, viscously driven 
accretion flow. Super-Keplerian rotation is also expected in 
the transition region between a standard Shakura & Sunyaev 
disc and an ADAF (Abramowicz, Igumenshchev & Lasota 
1997). We also note that the existence of super-Keplerian ro- 
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Figure 3. Distributions of angular momentum in the adiabatic models (7^ = 5/3) for four values of a (solid lines). The distribution of 
the Keplerian angular momentum is shown for comparison (dashed lines). For each a, the models have the values of the sonic radius Ts 
and the radial extension Tout = I'^out' given in Table 1. 



tation at the outer edge of some ADAF models can be seen 
in a few solutions previously obtained by Honma (1996). 

Just inside the outer boundary, the distribution of angu- 
lar momentum C is almost a power law with dln£/dlnr « 
0.77, larger than the Keplerian value of 0.5 (see Fig. 3). We 
note, that this value is practically independent of a and de- 
pends only weakly on r^, if rj^"^ ^ lOfg. In this case, the 
rotational energy of the accretion flow at the outer radii is 
being converted more efficiently into internal energy than 
in the case of a standard thin disc. As a result, both the 
temperature and the thickness of the disc increase rapidly 
inwards. A significantly sub-Keplerian rotation of this part 



of the disc implies that the radial pressure gradient becomes 
important in the balance of radial forces. 

Near to the black hole, solutions with high (super- 
Keplerian) and low (sub-Keplerian) angular momenta show 
qualitatively different behaviour. This was first noticed long 
ago by Abramowicz & Zurek (1981). When the bound- 
ary conditions are held fixed, high and low values of the 
angular momentum near to the black hole correspond to 
solutions with small and large values of a respectively. 
These two types of solution were referred to as 'disc-like' 
and 'Bondi-like', respectively. Differences between the disc- 
like and Bondi-like solutions in the special case of very 
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1 IQi 102 

r/r 

Figure 4. Distributions of radial velocity in the adiabatic models 

geometrically-thin accretion flows have been discussed by 
Muchotrzeb (1983), and later by several other authors, in- 
cluding Matsumoto et al. (1984), ACLS88 and Narayan et 
al. (1997). Our results show the existence of disc-like and 
Bondi-like solutions in agreement with the general results 
of Abramowicz & Zurek (1981). In particular, we confirm a 
more recent specific example of this general property, found 
by Narayan et al. (1997) in the special case of ADAFs; 
for a ^ 0.01, solutions are of the disc-like type, while for 
a 0.01 they are Bondi-like. For the disc-like solutions with 
small a, the sonic radius is located near to the marginally- 
bound radius, 2rg < rs ^ 2.3rg, and the distribution of 
angular momentum has a super-Keplerian part around the 
marginally-stable radius (see Figs 1 and 3). For the Bondi- 
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with 7g =5/3 (see the caption to Fig. 3 for further explanation). 



like type, corresponding to large a, the sonic radius is lo- 
cated closer to Vms, and the angular momentum is substan- 
tially sub-Keplerian everywhere in the region. Note, that our 
Bondi-like solutions have Vs < Vms- We did not find solutions 
with Ts > rms in case of a ^ 0.1. The existence of solutions 
with Ts ^ Vma IS questionable because the maximum ra- 
dius of the disc r^"^ decreases quickly with increasing sonic 
radius for any a (see Fig. 1). 

Differences between these two classes of solution can also be 
seen in the radial distribution of the pressure. Disc- like mod- 
els have pressure maxima outside the transonic point while 
these maxima disappear in Bondi-like models (see Fig. 5). 

Models with adiabatic index "fg — 4/3 only have prop- 
erties which are qualitatively similar to those with 7g = 5/3 
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Figure 5. Distributions of pressure in tlie adiabatic models with 7g =5/3 (see the caption to Fig. 3 for further explanation). 



if the viscosity is high [a J> 0.01). In particular, all of our 
high viscosity models with H/r < 1 and 7g = 4/3 are lim- 
ited in size. It could be, however, that there are also models 
which are unlimited in size but which our numerical proce- 
dure was unable to find: for purely technical reasons, it does 
not converge well for highly viscous models with large rout- 

Models with a ^ 0.001 show a very different behaviour, 
however. Models with a particular value of the sonic radius, 
rs = {rs)o are not limited in size and coincide asymptotically 
(for large radii) with the corresponding self-similar solutions 
of Narayan & Yi (1995). Models with > (rs)o are limited 
in size with r^"^ being a decreasing function of rs . For rs < 
(rs)o the solutions have H/r > 1 in the outer part and so 
do not satisfy the condition for the slim disc approximation. 



3.2 Bremsstrahlung models 

We now describe models with bremsstrahlung cooling, in 
particular, ones in which the inner ADAFs are matched with 
outer SLE type solutions. In the original SLE paper, ions and 
electrons were allowed to have different temperatures, and 
the assumption of equality of the local electron heating and 
cooling rates was used. We use a one-temperature approx- 
imation for plasma but, nevertheless, our bremsstrahlung 
solutions in the outer part can be treated as being of SLE 
type in the sense that they are locally cooled, i.e. F'^ « F~ . 

Fig. 6 shows the radial dependence of the relative thick- 
ness H/r, the importance of the advective fiux {F^ — 
F~)/F"'-, the angular momentum C, the equatorial pres- 
sure p, the Mach number A4, as defined by equation (2.24), 
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Figure 6. Selected properties of tlie bremsstralilung models with 7g = 5/3, a. = 10""^, M/M^ii^ = 10~^ and M = WMq. The dotted, 
dotted-short-dashed, short-dashed and solid lines represent the cases Ts/rg = 2.09, 2.1, 2.2 and 2.7, respectively. The behaviour of the 
analytic SLE solution (3.1) is shown for comparison by a dotted-long-dashed line in the aspect ratio H/r plot. A dotted-long-dashed line 
in the specific angular momentum C/vgC plot shows the Keplerian angular momentum distribution. The bottom left frame shows the 
run of the Mach number M as defined by equation (2.24). The global solutions consist of two different parts with a relatively narrow 
transition region. At small r, the solutions behave like ADAFs, while at large r they behave like one-temperature SLE solutions. The 
model with Ts/vg = 2.7 (solid line) has special properties: it consists of the SLE solution everywhere except for the region inside 3rg. 
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and the velocity V for models with 7g = 5/3, a = 10 and 
M/MEdd = lO"'^. Here Msdd = ATrGMmp/car is the Ed- 
dington accretion rate. The black hole mass was taken to be 
M = IOMq. In these models we varied the position of the 
sonic radius. Dotted, dotted-short-dashed, short-dashed and 
solid lines represent the cases rs/vg = 2.09, 2.1, 2.2 and 2.7, 
respectively. Inward of the transition radius, the solutions 
become ADAFs with (F+ - ~ 1, and they have 

all of the properties of the adiabatic models discussed in the 
previous paragraph. The transition radius between the two 
types of solution depends on r^. This dependence, as well as 
the limitations of the minimum value of r^, are very similar 
to those for r^"^ (see Fig. 1, square dots) and {rs)crit in 
adiabatic models. The model with rs/vg = 2.7 (solid lines in 
Fig. 6) has special properties, because it becomes an ADAF 
just inside 3rg, and the transition radius is very close to the 
sonic radius. 

Outside the transition radius, each solution matches 
smoothly with an SLE type solution, which is unique for 
given Q, M and M, and does not depend on the position 
of the transonic point or the location of the transition 
radius. Note that we have not assumed any extra physical 
effect that triggers the transition. The transition occurs due 
to the standard accretion processes, described by the stan- 
dard equations for slim accretion discs. 

Of course, one my construct SLE solution analytically, tak- 
ing Q = flK, = F~ and zero viscous torque at r = 3rg, 

where rh = M/MEdd- The analytic SLE solution (3.1) is 
shown in Fig. 6 by the dotted-long-dashed line. There is 
good agreement between the analytic solution (3.1) and the 
numerical solution outside the transition radius. Both solu- 
tions have similar slopes inside ~ W^rg. Fig. 6 shows that 
the SLE type solutions have negative [F'^ — F~)/F^ (see 
the discussion by Abramowicz 1996). In the standard SLE 
solution F'^ — F~ . 

In Fig. 7 we show properties of the 7g = 4/3 models. 
Short-dashed, dotted, long-dashed, dotted-short-dashed and 
solid lines represent the cases rs/rg = 2.13075, 2.131, 2.14, 
2.255 and 2.7, respectively. All of the models experience the 
ADAF-SLE transition. Note, that in models corresponding 
to the short-dashed and dotted lines, the ADAF part is very 
extended. In contrast to the previous (7^ = 5/3) case, we 
found an unphysical behaviour of the solutions for r > 10^ rg. 
In this radial range, the solutions have the angular momen- 
tum jC decreasing outwards, which is dynamically unstable. 
They also deviate strongly from the thermal structure char- 
acteristic of the SLE solution, having the cooling rate F~ 
significantly larger than the heating rate F"*". This behaviour 
does not depend either on the choice of the location of the 
outer numerical boundary Vout, or on the value taken for 

^out • 

4 DISCUSSION AND CONCLUSIONS 

We have constructed numerical solutions for two classes of 
models of optically thin accretion discs around black holes: 



adiabatic models, with no radiative cooling, and models with 
bremsstrahlung cooling. 

All of the adiabatic models are of the ADAF type. Those 
with 7g = 5/3 and H/r ^ 1 everywhere are limited in size: 
their radial extensions are smaller than about 10^ rg. Inde- 
pendently of whether viscosity is high or low, the vertical 
thickness H of these models tends to zero at the outer lim- 
iting radius. Models with 7g = 4/3 and H/r ^ 1 extend to 
radial infinity only when the viscosity is low (we constructed 
models with a ^ 10~^) and when the location of the sonic 
point has exactly a particular value Vs — (t's)o. Obviously, 
this is the eigenvalue of the problem. It is therefore not sur- 
prising that these models have, at large distances, the same 
properties as the models corresponding to the self-similar 
solution of Narayan & Yi (1995). For > (rs)o our models 
are radially limited. For discs with high viscosity (a <: 10"^) 
we could not find disc-like solutions with unlimited radial 
size having H/r ^ 1 everywhere. We cannot conclude that 
such solutions do not exist, due to numerical problems which 
our code experiences in this case. Radially unlimited solu- 
tions that we found in this range of parameters, all have 
the maximum value of H/r being an increasing function of 
the location taken for the outer numerical boundary rout- 
In particular, for large rout one always has H/r S> 1. This 
indicates that the one-dimensional approach based on verti- 
cal integration which assumes H/r ^ 1 is not self-consistent 
and slim models for this type of disc cannot be trusted. 

Bremsstrahlung models consist of two parts which are 
joined by a smooth transition region. In some cases the tran- 
sition region is relatively narrow (models with 7g = 5/3), but 
in other cases it can be larger than the ADAF. Inside the 
transition radius, the solution is very similar to the ADAF- 
type adiabatic solution. Outside it, the solution is of the SLE 
type, with the entropy of the matter increasing with radius. 
Our numerical SLE-type solutions coincide closely with the 
analytical ones. Obviously, though, this kind of global solu- 
tion cannot exist in nature, because the outer SLE-type disc 
is thermally unstable. 

The most interesting property of the ADAF solutions 
found in this paper is their behaviour near to the transition 
radius. It was known previously (see e.g. discussion given by 
Chen et al. 1995) that the ADAF-type accretion flows can- 
not extend to arbitrarily large radii — they may extend out 
to ~ lO^Tg. The limitation is due to the fact that at large 
radii the accretion time scale exceeds the free-free cooling 
time. This case is represented by the model shown by the 
short-dashed line in Fig. 7. At the range of radii from lO'^rg 
to around lO^rg the short-dashed model is an ADAF, which 
is similar to the self-similar solution of Narayan & Yi (1995), 
where H/R ~ const and Q/i^K ~ const- At r ^ lO^^g the 
radiative cooling becomes important in the energy balance 
and the solution ceases to be of the ADAF type. However, 
the end of the ADAF-type part for the other models shown 
in Figs 6 and 7 must be of a different nature. It cannot be 
connected to the cooling mechanisms, because both adia- 
batic and non-adiabatic models show very similar behaviour. 
One should note that the position of the end of the ADAF 
is determined by the value of the parameter rs'- a larger 
corresponds a smaller transition radius. 

The absence of the radially-unlimited ADAF solutions 
for 7g = 5/3 could be connected with the fact that there are 
no disc-like asymptotic self-similar solutions for 7g = 5/3. 
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The case jg = 5/3 is singular (see Narayan & Yi 1995). 
On the other hand, we note, that the analogous self-similar 
solutions derived from the slim disc equations in the form 
given by Chen et al. (1997) are free from the jg — 5/3 
singularity. Our slim disc equations are similar to the ones 



used by Narayan & Yi (1995), and therefore it is natural 
that our results agree with theirs. The question, however, 
remains of why different versions of the slim disc equations 
produce qualitatively different solutions at large radial dis- 
tances r ^ lO^r-g? Are seemingly 'formal' details of the ver- 
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tical integration crucial for the physics of slim disks? Is the 
slim disc approach questionable in application to ADAFs? 
The best way to answer these questions would be by con- 
structing fully two dimensional, non-stationary models of 
ADAFs. 

There is a less important, but related, problem here. We 
have not obtained global solutions with shocks like those re- 
ported by Chakrabarti and collaborators (see Chakrabarti 
1996) despite the very similar physics and boundary condi- 
tions assumed. Some of our solutions have transition regions 
which are narrow, but smooth with their inner structure be- 
ing numerically well resolved. These narrow transitions are 
definitely not shocks. It is likely that this will be taken as 
an argument in the recent discussion 'shocks or no shocks' 
but we do not want to join in with this ourselves because 
we are convinced that the discussion is artificially inflated 
and of only academic interest. It does not touch a really 
interesting point here that one should be aware of. Several 
authors have recently studied general properties of transi- 
tion regions between ADAFs and thin discs (see e.g. Regev, 
Lasota & Abramowicz 1997; Abramowicz et al. 1997b) and 
concluded that the physical conditions in the transition re- 
gion are very complex, and cannot be described by any sim- 
ple one-dimensional, stationary approach. Details of the ver- 
tical structure and dissipative processes are fundamental for 
the structure of the transition region and must be fully taken 
into account. Thus, one should first construct a more physi- 
cal model of the transition region, before speculating about 
its nature. 
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